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We present a realistic dark energy model derived from particle physics. Our model has essentially 
no free parameters and has an equivalent fit to the observational data (CMB, SNla and LSS) as 
LCDM and a better fit than the best effective wi^z) model. With the lack of a clear determination 
of the cosmological parameters theoretical considerations should be taken seriously to distinguish 
between dark energy models. 
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Recent cosmological observations imply the uni- 

verse is dominated by a dark energy fluid. How- 
ever, there is no theoretical understanding what this 
dark energy is. The simplest and most common as- 
sumption is a cosmological constant model. Compar- 
ing LCDM to wmap, SNla and LSS data one finds 
a large degeneracy among the cosmological parameters 
Pq = {Ho,i^DE,^bh^ ,T,ns, As) I3|-@- Several attempts 
have been made to infer from the data the best dark en- 
ergy model by parameterizing the equation of state pa- 
rameter w{z) as a function of rcdshift z 5j. The analysis 
shows that the degeneracy of the cosmological parame- 
ters Pq remains but the central values varies. The best fit 
model for a dynamical 'w{z) and LCDM have an equiva- 
lent x^/dof = 1.08 value |4[,[| but LCDM has a better 
Akaike or Bayesian information criteria [^. With the 
lack of a clear determination of Pg theoretical consider- 
ations should be taken seriously to distinguish between 
dark energy models. 

Even though the LCDM model has a good fit to the 
cosmological data it is hard to accept it from a theoreti- 
cal point of view. It does not answer some fundamental 
questions as where it comes from and why does it have at 
present time an energy of the same order of magnitude as 
dark matter, i.e. the coincidence problem. From a theo- 
retical point of view we would like to have a dark energy 
model which is well motivated from particle physics, an- 
swers the fundamental question which a LCDM does not 
and has at least an equivalent fit to the data. The best 
theoretical model should have the least number of free 
parameters and should give reasonable explanations and 
values for them. 

In this work we will present a realistic dark energy 
model derived from particle physics. We will show that 
this model has, indeed, an equivalent fit to the observa- 
tional data as LCDM. It has essentially no free param- 
eters so the cosmological Pq quantities are the only free 
parameters. For most of the cosmological time the energy 
density is proportional to radiation since the particles are 
massless. It is the dynamics of the model that give rise to 
a late time phase transition at a scale Ac = 0(10— lOOey) 
and an effective scalar field appears which gives the dark 
energy today. For energies larger than Ac the particles 
are massless and rcdshift as radiation and below Ac the 
particles acquire a non-perturbative mass. The model is 
based on gauge theory and has an equivalent structure as 
the well established standard model of particle physics. 



The dark energy gauge group gives an effective poten- 
tial of the form V oc 4>~", where 4> is the dark energy 
(quintessence field) and the power n is fixed by the num- 
ber of fields. Inverse power low potentials "IPL" were 
first introduced by and ^ to use ADS superpoten- 
tials to get an effective IPL potential. However, the 
models of w'l had no particle physics interpretation and 
the models of ^^J had n > 2 and are phenomenological not 
viable. We were the first ones to study phenomenologi- 
cal acceptable ADS superpotentials giving an IPL with 
n < 2 and our model has n ~ 2/3 JJ]]. We would like to 
emphasis the following theoretical points of our model: 

• The potential is calculated from our best theoretical 
theories, i.e. gauge theories 

• The superpotential is exact and the potential is stable 
against quantum corrections 

• The appearance of quintessence field is due to a late 
time phase transition given by Ac 

• The condensation scale Ac is determined by first prin- 
ciples (it is given in terms of the number of fields) 

• The initial energy is given also only by the number of 
fields 

• The number of fields is determined by imposing gauge 
coupling unification 

• The solution is an attractor even though the 
quintessence field has not reached it yet 

• Our dark energy model has no free parameters 

• It has an equivalent fit to the observational data (CMB 
and SNla) as the best fit model 

• It "solves" the coincidences problem 

This last point is due to a late time phase transition. 
The onset of the quintessence field is at about matter- 
radiation equality since the scale factor Oc (i.e. at Ac) 
is in cosmological times very close to matter-radiation 
equality Ueq ( ac/aeq — 10~^). This means that poc/Pr 
is constant from the end of inflation ainf — 10~'^° until 
flc — 10^^ which accounts for most of the time. 

The dark energy model is simply a SU {Nc = 3) gauge 
group with Nf = 6 elementary particles in the funda- 
mental representation and with only gravitational in- 
teraction with the standard model of particle physics. 
The phase transition scale Ac is determined only by A^c 
and Nf and the value of the gauge coupling constant at 
some arbitrary energy scale. Motivated by string the- 
ory our dark energy model is constrained to be unified 
with the standard model gauge groups at the unifica- 
tion scale. For viable cosmological models gauge con- 
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pling unification determines the values of Nc,Nf [To| . 
With this non-trivial constrained the evolution of the 
gauge coupling constant and Ac are completely fixed. At 
high energies the dark elementary fields are massless and 
Pdg oc a~^, pdg/Pt is constant and the ratio is given 
only in terms of the number of particles (c.f. eq.J^J). 
At lower energies a phase transition takes place due to 
a strong gauge coupling constant. At this scale the dark 
elementary fields are bind together producing gauge in- 
variant states. The relevant scale for this process is the 
condensation scale Ac and for a gauge group SU{Nc) 
with Nf matter fields in susy is given by the one-loop 

renormalization group equation Ac = ^gut&~ "Sg^t 
where bo = 3A^c — Nf is the one-loop beta function and 
Agut — lO^^GeV, Qgut — \/47r/25.7 are the unification 
energy scale and coupling constant, respectively. In our 
dark energy model we have Nc = 3, Nf = 6 giving bo ~ 3 
and Ac = A2eV. Strong gauge interactions produce a 
non-perturbative potential V. This potential can be cal- 
culated using ADS potential The superpotential for 
a non-abelian SU{Nc) gauge group with Nf massless 
fields isW = {N^ - Nf){Kl- ldet{QQ)Y/^^''-^f\ The 
scalar potential in susy for one dynamical meson field 
0, the pseudo-Goldstone boson, is V — e'^ ^^IW^P with 
= dW/d(l), giving |3,[l3 



(1) 



with c = 2Nf = 12. We have included the exponential 
factor in eq. because once gravity is taken into account 
it gives the correct normalization in the Einstein frame. 
For values of the field much smaller than the Planck mass 
rup = 2.4 10^* GeV = 1 the exponent term is irrelevant 
but for ~ 1 it gives a (small) correction. The evolution 
is basically determined by the inverse power law expo- 
nent . In order to study the evolution of eq.QJ the 
initial condition on (j) must be set and the only natural 
initial value is 4>i = Ac since it is precisely Ac the relevant 
scale of the physical binding process. The potential in 
eq. has a minimum aX (f) = ^/n = ^2/3 and the value 

of the potential is F|™„ = c2Ac"^^/^ei/3(3/2)i/3. The 
evolution of with different initial conditions is shown 
in fig|21 Even though the evolution of (p depends on Jl^; 
its value at present time Oo = 1 is (almost) independent 
of the initial conditions. Global symmetries and susy 
protect the mass of the quintessence field (j). In fact the 
ADS superpotential is exact and receives no corrections 
. The only radiative corrections arise due to the Kahler 
potential. However, our model is stable against radiative 
corrections the shape of the potential remains the 
same and the v.e.v. of suffers a small shift but this does 
not affect the cosmological predictions since the solution 
is an attractor. The complete evolution of the dark en- 
ergy group can be seen in figH Notice that pdg tracks 
radiation for a long period of time, including nucleosyn- 
thesis "NS" epoch, since all the particles are massless 
and the onset of the quintessence field is at a very late 
time (oc — 10~^) and close to matter-radiation equality 
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FIG. 1: We show the evolution of pr,pm and poG (dotted, 
dashed and solid lines, respectively) as a function of Logio{a) 
with = 1 with no extra degrees of freedom, i.e. ~ 0. The 
value of equivalence is n^q = — 3.6 while the phase transition 
takes place at = —5.9. 
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FIG. 2; We show the evolution of for diflerent initial 
conditions floGi = 0.01,0.03,0.06,0.11 (dotted, dashed, dot- 
dashed and solid lines, respectively). 



{tteq — 10""*). The fact that Ac is so small and the ap- 
pearance of the quintessence field is at such a late time 
"solves" the coincidence problem since it implies that an 
accelerating universe will necessarily be at a scale factor 
larger than Oc and Oeg, close to present time. In our case 
the scalar field is present only at a > Oc and it does not 
track the background at all, as can be seen from fie:s lll3l 
and El 

The total energy density of the universe is given px = 
Psm + Pdg + Ph, where psm is the energy density of the 
standard model (we will take it as the minimal super- 
symmetric standard model) and we have allowed for ph 
hidden particles which could be present at earlier times. 
Examples of these particles are the ones requires to break 
susy in the standard model and are widely predicted by 
string theory. The only constrained we imposed to ph 
is that they do not contribute to the present day energy 
density nor at NS. One important constrained on the en- 
ergy density at earlier times is given at NS. The extra 
amount of relativistic energy density at NS is usually pa- 
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rameterized by an effective neutrino AA''^ < 0.2 — 0.3 
which impUes an upper bound fl^xiNS) < 0.045 in the 
most stringent case and flexiNS) < 0.09 at two a [T^ . 
However, NS is still fine with AA'^i^ < 4.1 which implies a 
bound riexiNS) < 0.4 if new physics is involved at that 
time (as neutrino asymmetry) |l3| . 

Our dark energy model has no free parameter, how- 
ever, the total energy density could in principle con- 
tained more particles than the standard model and our 
dark group, i.e. gn 0- In fact extra particles would 
be needed to break susy in the visible sector. To sim- 
plify the discussion we will assume that these particles, 
if present, acquire a large mass (larger than TeV) so that 
they do not affect NS and that they are coupled to the 
standard model at high energies. Examples of this are 
the gauge mediated susy breaking models, extensively 
studied in particle physics context, which require a en- 
ergy scale Ah = O{10^~^)GeV in order to give masses 
of the order of TeV to the standard model particles ■ 
Since the hidden sector is only coupled to our dark model 
via gravity, it will induce a susy breaking to the dark 
group of the order nig = A^/m^ which is smaller than Ac 
and does, therefore, not affect the running of the gauge 
coupling and the onset of our potential given in eq.l^. 
These extra particles will not contribute to the energy 
density at NS nor at present time but they are impor- 
tant for determining the ratio of temperature between 
the standard model and our dark group. The effect of 
these extra particles is to decrease the ratio of the tem- 
perature and ^DG- From entropy conservations one can 
determine ^lDG giving 11] 
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Td 
T 



gdec 



1/3 



gsrn +gQ{TD/Ty 

(2) 

with q = 10.75, gs,n = 10.75 at NS and g = 2 x 10.75 x 
4/11, gsrn = 3.36 at Ac and gdec = gnsSM + gh, with 
gMSSM — 228.75 and gh the degrees of freedom of the 
MSSM and the hidden sector, respectively. The term 
4/11 in q at Ac takes into account for neutrino decoupling 
and 10.75, 3.36 are the relativistic degrees of freedom of 
the standard model at NS and at present time while gq — 
97.5 are the degrees of freedom of our model ITT]. An 
upper bound on dq , can be established if there are 
no hidden particles, gt = 0, and eq.© gives 



VLDoiNS) < 0.13 



^^DG^(Ac) < 0.11 (3) 



where i^DGi stands at Ac, i.e. the onset of the 
quintessence field. Notice that the upper limit is still 
within the existing NS bounds (not the most stringent 
one) since AiV^ < 1 implies ilex < 0.14. For g^ 
than we get a smaller V,j:)g{NS) and for Qdg{NS) — 
0.09, 0.045 we require gh = 90, 327, respectively, giv- 
ing i^DGi — 0.076, 0.037. Typical gauge mediated susy 
breaking models have gh — 0(200) [iJl- Clearly we have 
no precise knowledge of how many hidden sector particles 
are nor how susy is broken. Furthermore, the contribu- 
tion of gh would be the same for all dark energy models 
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FIG. 3: We show the evolution of w as a function of n = 
logio{a) from Uc (i.e. at Ac) to present day no — for differ- 
ent initial conditions QoGi = 0.01, 0.03, 0.06, 0.11 (dotted, 
dashed, dot-dashed and solid lines, respectively). 
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FIG. 4: We show the late time evolution of lii as a func- 
tion of redshift z for different initial conditions QdGi = 
0.01, 0.03, 0.06, 0.11 (dotted, dashed, dot-dashed and solid 
lines, respectively). 



and the uncertainty lies in the hidden sector and not on 
our dark group model which, as stated above, has no free 
parameter. 

In fig|3|we show the evolution of w below Ac, i.e. for 
a > flc (above Ac it is simply w = 1/3). It has an initial 
kinetic period w = 1 following a region with w — —1, 
lasting almost the same amount of e-folds as the first one, 
and it grows later to its present value Wo- We see that 
the value of Wo depends slightly on the initial conditions 
^DGi- At equivalence LogiQ[aeq\ — —3.6 the value of 
dark energy is ^Idg — 10"^ while at last scattering (a ~ 
1/1090) it is ^DG - 10"'^, quite smaU in both cases. In 
fig 01 we show how w increases as a function of redshift z 
for different values of QoGi- In figEt we plot the present 
day values of Wo and w'^ = dw/dz\o for different initial 
^DGi and final QoGo- Notice that the variation is very 
small. In fact we have —l<Wo< —0.92 and —0.06 < 
< 0.19 if ^DGo varies from 0.65 < iloGo < 0.75 and 
^DGi < 0.11. For larger flDGi > 0.1 we get a flat and 
small Wo, very close to -1, with a slightly negative w'^ 
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FIG. 5: We show the values of w and its derivative w' — 
dw/dz\o at present time for different values of Qqo and 
f^DGi- The dotted lines have constant Qqo and from bot- 
tom to top they have Q.qo = 0.66,0.68,0.70,0.72,0.74, re- 
spectively. The solid lines have QoGi constant with Qoa ~ 
0.0002,0.01,0.02,0.03,0.04,0.05,0.06 from top to bottom. 
Notice the resulting w, w' region is very small and it has an 
upper limit to w and w' . 
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TABLE I: We show the best fit model compared to the 
Golden SNla set for different values of Q,DGi and LCDM. 



while for values ^oa < 0-1 we have a small positive 
w'^ and a larger Wo- As seen from fig|SJi the variation is 
more sensitive with respect to iloGi than to floGo- The 
values of Wo, w'^ are insensitive to Ho- At this point, we 
would like to emphasis, again, that ^dgi is not a free 
parameter, we do not marginalize it, but we show the 
slight dependence of Wo, w'^ over it. 

We now compare our model with the data. The 
likelihood contours for the Golden SNla set of 
^DGi vs Hdgo are shown in figEI The best fit to the 
Golden SNla set for different initial conditions are shown 
in tabled We see that for the different values of fi^i we 
get an equivalent fit as a LCDM. The acceleration red- 
shift (where the universe begins to accelerate) is centered 
around Zac = 0.65 ± 0.05 with xLia < 178. The best fit 
model has a very flat w with —1 < Wo < —0.94, even 
though the ^Idgi = 0.01 case with Wo = —9.94, w'^ — 
0.15 does a reasonable job. 

To compare our model to the wmap data we use KINK- 
FAST [ll (a modification of CMBFAST [l6| to include 
a dynamical dark energy). The best fit to wmap of our 
dark energy model and LCDM can be seen from table 
im The parameters , t are degenerated and different 
values would give the same Xtot 
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FIG. 6: We show the likelihood contours as a function of 
Q,DGi,^Qo compared to the Golden set of SNla. We have 
xLia < 178, 178 < xLia < 179, 179 < xLia < 180 and 
180 < Xsnia < 181 (light blue, red, light green and black, 
respectively) . 
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TABLE II: We show the best fit model compared to wmap 
data for different values of Qdgi and LCDM. 



Combining the Golden SNla set with wmap we show 
in tablc lllll the best fit model for different values of floGi 
and for the best LCDM. Notice that all four initial condi- 
tions have an equivalent fit and therefore, the data is not 
sensitive enough to distinguish between them. The total 
number of degrees of freedom is 157 for the Golden SNla 
set and 1342 for the wmap data. We get for our best fit 
model a xL = 1605.8 with xLt/dof = 1.08. The dof 
is the same for our dark energy model as for LCDM. It 
is interesting to note that after an exhaustive analysis of 
model independent evolution of dark energy ^ the best 
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TABLE III: We show the best fit model compared to the 
Golden SNla set and wmap for different values of QoGi and 
a LCDM. All models have Us = 0.97, wb = 0.0227. The 
number of degrees of freedom is the same for our dark energy 
model and LCDM, with a total dof 13424-157-6= 1493, giving 
a xlt/dof = 1.08. 
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fit model has Xtot — 1602.9 at the price of introducing 4 
extra parameters giving the same X^ot/'^'^f — 1-08 but a 
worse Akaika or Bayesian criteria Q ■ The Akaika infor- 
mation criteria ^3 requires the smallest AIC = 'x^^^ + 2k 
while the Bayesion [13 has BIC — Xtot + ^ InN, where 
k is the number of paramters and N the number of 
data points. Our dark energy model and LCDM have 
AIC = 1605.8 + 12 = 1617.8, BIC = 1664.3 while the 
best w{z) has AIC = 1622.9, BIC = 1676. We see that 
for both criteria our dark energy models has a better fit 
than the best effective w{z) model and an equivalent fit 
as LCDM. The evolution of our dark energy model and 
CTg lie at the central values obtained in 

As a matter of completeness we have determined the 
best fit if the exponential factor in eq.(Q) is not present 



(i.e. for V = c^A 



(69, 0.71, 0.0228, 0.12, 0.97) giving 
= 177.7 + 1428.9 ^ 1606.6 with 



nDG^ - 0.11 a P5 

^tot Xsnla ~^ ^wr 

Xtot/^of — 1-08 ^"^d an equivalent fit as before. 

To see the difference between our dark energy model 
compared to LCDM we have calculated the CMB us- 
ing the same Pg in both cases. For small Hdgi (i-e. 
^DGi = 0.01, 0.03) we get a variation of about 0.2 — 0.5% 
increase in the three peaks while an increase of about 
1 — 3% at lower multipoles with a resulting Ax^^ap — 10 
in favor of our dark energy model. On the other hand 
for i^DGi — 0.06 or 0.11 the variation is much smaller, 
lees than 0.05%, giving, therefore, an equivalent Xwmap- 
However, by changing Pg we can get an equivalent fit 
for LCDM as for our dark energy models. The ISW ef- 
fect for a dynamical w(z) gives a nontrivial (but small) 
imprint on the power spectrum ^ but this contribution 
can be taking into account by a change in Pg . The great 
degeneracy on the cosmological parameters Pg reduces 



the possibility of distinguishing between different dark 
energy models at present time. The LSS data do not 
place a stronger constraint on the evolution of dark en- 
ergy than the SNla and wmap data 0-0- Without an 
independent way of measuring Pg it is hard to distinguish 
from the data between different dark energy models and 
we should choose the best theoretical model. 

Let us conclude. We have presented a realistic dark en- 
ergy model derived from particle physics that has essen- 
tially no free parameters. The reason why the universe 
accelerates at such a late time is because the dark en- 
ergy (quintessence field) appears at a very late time due 
to the gauge group dynamics. Our model has a almost 
flat w, with — 1 < Wo < —0.94, and has a better fit to 
the cosmological observations than the best w{z) effective 
model and an equivalent fit than LCDM. Unfortunately 
the great degeneracy on the cosmological parameters Pg 
reduces the possibility of distinguishing between different 
dark energy models from the observations. Therefore, we 
are left, at present time, with the most appealing theoret- 
ical model as the best candidate for dark energy. This is 
the model with the least number of free parameters and 
with the best explanation of the origin of dark energy 
and its late time appearance. Our dark energy model 
satisfies these criteria. 
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